Temperature and Disorder Chaos in three-dimensional Ising Spin Glasses 
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We study the effects of small temperature as well as disorder perturbations on the equilibrium 
state of three-dimensional Ising spin glasses via an alternate scaling ansatz. By using Monte Carlo 
simulations, we show that temperature and disorder perturbations yield chaotic changes in the 
equilibrium state and that temperature chaos is considerably harder to observe than disorder chaos. 
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The fragility of the equilibrium state of random frus- 
trated systems such, as the Edwards- Anderson Isiiig spin 
glass has been predicted a long time ago Jl,[5| and 

analyzed on the basis of scaling arguments [&, '7] . These 
scaling arguments predict that the configurations which 
dominate the partition function change drastically and 
randomly when the temperature or the disorder in the in- 
teractions between the spins are modified ever so slightly. 
The temperature chaos and disorder chaos effects have at- 
tracted considerable attention both from theory and ex- 
periment because of their potential relevance in explain- 
ing the spectacular rejuvenation and memory effects ob- 
served in hysteresis experiments in spin-glasses @, 0, 
as well as other materials, such as random polymers and 
pinned elastic manifolds. Although there is evidence of 
disorder chaos in spin glasses, temperature chaos remains 
a controversial issue ^Jj, il^, dS,, lM, dM ' whereas for ran- 
dom polymers or pinned elastic objects [l6| chaos in gen- 
eral is well established 



model and Monte Carlo methods used, followed by the 
disorder- and temperature-chaos scaling approaches. We 
conclude with the results of our simulations of the 3D 
Ising spin glass and a general discussion. 

Model and numerical method The Edwards- Anderson 
1] Ising spin glass is given by the Hamiltonian 



n 



ijSiSj, 



(1) 



where the Ising spins Si G {±1} are on a cubic lattice 
with N = vertices and the interactions Jij are Gaus- 
sian distributed random numbers with zero mean and 
standard deviation unity. The sum is over nearest neigh- 
bor pairs. The model undergoes a spin-glass transition 
at Tc = 0.951(9) 0, [13, 111- 

The order parameter of the system is defined via the 
overlap between two copies a and i.e., 
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Following previous studies [l^ . [T^ . we probe tempera- 
ture chaos when the temperature between both replicas 
is shifted by an amount AT. Disorder chaos is studied 
by introducing a perturbation AJ in the disorder, i.e., 



Despite this lack of consensus, it has been surmised 
that temperature chaos would only be observable in spin 
glasses at very large system sizes and for large temper- 
ature changes 2^, 21 1 thus making its presence unfath- 
omable in simulations. These claims have recently been 
challenged. In particular, recent results point towards 
the existence of temperature chaos in four-dimensional 
Ising spin glasses 22] where the free energy of a domain 
wall induced by a change in boundary conditions changes 
its sign chaotically with temperature in accordance with 
the droplet/scaling theories 0, @, In this work we 
study the overlap between states at different tempera- 
tures and disorder distributions directly for a physically 
relevant three-dimensional (3D) Ising spin glass. Our 
results show that the scaling laws that arise from the 
droplet theory are indeed well satisfied in 3D provided 
low enough temperatures are considered, although small 
corrections need to be applied. In addition, we show that 
temperature and disorder chaos have similar scaling func- 
tions. By rescaling the characteristic length scale in the 
problem, we show that disorder chaos appears at much 
shorter scales than temperature chaos 

The paper is organized as follows: We discuss first the for disorder chaos, respectively. 
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which leaves the disorder distribution invariant. In 
Eq. Q x is a Gaussian distributed random number with 
zero mean and standard deviation unity. To monitor the 
changes induced by the perturbations of the system, we 
compute the chaoticity parameter 1^ given by 



Qti,i 



. \J ^1Ti,Ti){1t2,T2) 

for temperature chaos and by 



(4) 
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TABLE I; Parameters of the simulation for each system size 



is the number of samples and N^^ is the total 



number of Monte Carlo sweeps used for equilibration for each 
of the 2Nt parallel tempering replicas for a single sample. 
An equal number of sweeps is used for measurement. The 
minimum temperature simulated is Tmin = 0.20, the highest 
Tmax = 2.0. For disorder chaos the disorder shifts used in 
Eq. © are AJ = 0.001, 0.005, 0.02, 0.05, 0.1, 0.2, and 0.5. 
For temperature chaos we compute the overlap between Tmin 
and Ti with i £ {2, . . . , Nt} 29]. 
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In Eqs. (j4|) and ([5]) is the square overlap, Eq. be- 
tween two copies at different temperature/disorder. Here 
(• • • ) represents a thermal average and [• • • ]av represents 
a disorder average. In order to access low temperatures 
necessary to probe temperature chaos, we have used the 
parallel tempering 2^, 27 1 Monte Carlo method in com- 
bination with the equilibration test presented in Ref. [28| . 
Simulation parameters are listed in Table [D 



Disorder and Temperature Chaos In what follows we 
discuss how chaos can arise in spin glasses using the early 



[6L 0], and [3d]. Within 
|3l|, the low- lying exci- 



arguments presented in Refs. [4 L 
the droplet theory framework 6 
tations above the equilibrium state are obtained by flip- 
ping compact connected clusters of spins called droplets. 
A droplet of size ^ has a fractal surface of dimension 
ds < d, and its excitation free energy > is distributed 
via PTiF,e) = ['y(T)£]-'^p[F/-f{T)e'^], where p{x) is a 
scaling function assumed to be nonzero at a; = and 
which decays to zero for large x. The free-energy ex- 
ponent is argued on general grounds to be such that 
< 6 < ds/2 and 7(T) is the free-energy stiffness (which 
goes to zero at Tc). The droplet's entropy S can be writ- 
ten as ^ = cr(T)£'^=/2, where a is the entropy stiffness. 
Temperature chaos appears if the free energy of a droplet 
changes its sign when the temperature is modified. As 
noted in Refs. Q and the length scale at which this 
happens can be estimated by noting that the energy of a 
droplet does not change much with temperature. There- 
fore, if one considers a droplet at temperature Ti with 
free energy F{Ti), then at temperature T2 > Ti 



F{T2) « F(Ti) + ri5(Ti) - T2S{T2) 



(6) 



Because for typical droplets F{Ti) = 'yiTi)f and S = 
cr(T)£''''/^, the free energy excitation of such droplets be- 
comes generally negative at temperature T2 (so that the 
droplet has to be flipped) for length scales larger than 



the chaotic length @, 0] defined as 
7(Ti) 



4 = 



T2 <J{T2)-Ti a{Ti) 



1/C 



with C = — — 
^ 2 



(7) 



Usually, small temperature changes are studied such that 
S{Ti) « S{T2). Here however, we do not use this approx- 
imation and since in the low temperature phase, when 
one can define droplets, the entropy is proportional to 
Slii, we write 



3/2 



-1/C 



(8) 



Equation ([5]) shows that when temperature is changed, 
equilibrium configurations are changed on scales greater 
than £c- Notice that by keeping the temperature depen- 
dence of the entropy, we obtain a slightly different scaling 

than usually considered [6, 7], where a factor AT appears 

3/2 3/2 

instead of r2 — . While this makes no difference 
for small AT (which is the case in all simulations per- 
formed so far), this can be significant for temperature 
differences larger than the ones considered in this work. 
Similar arguments can also be applied to the case of a 
random perturbation in the disorder @, 0] ~ see for in- 
stance Ref. [33] - where one obtains 



4(A^) oc AJ" 



-1/C 



(9) 



Considering system-size excitations, these arguments 
thus suggest that the chaoticity parameters defined in 
Eqs. (IH) and (O have the following scaling behavior 



QTi,T.(L,Ti,T2) = FT[i/4(Ti,r2)], 
Qaj{L,AJ) = Fj[L/4(AJ)], 



(10) 



where F{x) is a function with F{0) = 1 that decays at 
large x. In what follows we test the aforementioned scal- 
ing relations via Monte Carlo simulations. 

Numerical results The behavior of the chaoticity pa- 
rameter for disorder chaos at zero temperature has re- 
cently been studied in Ref. (sl] (see also Ref. [11]) and 
perfectly satisfies the predictions of the droplet model. 
We thus concentrate here on finite temperatures where 
the scalirig relations were only tested in the simulations 
of Refs. [l3|, and [S^l in two and four space di- 

mensions. We have performed low-temperature Monte 
Carlo simulations of the 3D Ising spin glass (see Table 
m. As can be observed in Fig. (TJ scahng the data for 
disorder and temperature chaos according to Eqs. (fTU]) 
works extremely well. The best scaling collapse deter- 
mined by a nonlinear minimization routine [25| yields 
C « 1.04 for temperature chaos and ( « 1.16 for disor- 
der chaos, which is in rather good agreement with the 
accepted value C ~ 1-1 from 4 « 2.6 [21, Hi] and 9 w 0.2 
31, m, [111; see Eq. Q. Notice that we have a good 



scaling of the data even when T2 is larger than Tc. We 



3 




temperature 
f = 1.04(10) 

5 10 15 20 Zb 



m 



0.8 



L 

□ 4 
O 5 
O 6 - 
A 8 
010 



disorder 

f = 1.16(10) t 

. .ill. I. I. I. il. mil. 

12 3 4 



5 10 15 20 25 1 2 3 4 5 

FIG. 1: (Color online) Scaling plot of the chaoticity parameter 
for temperature chaos (left panel, T\ = 0.20) and disorder 
chaos (right panel, T — 0.2) using L/£c as a scaling variable. 




L/iK) i^/iK) 
FIG. 2: (Color online) Scaling plot of both disorder and tem- 
perature chaos. Left: Chaoticity parameter Q as a func- 
tion of L/£c. Scaling performed with ( = 1.16. Note that 
£c{AT) ~ A£c{AJ) with A 8.7, i.e., disorder chaos appears 
at considerably shorter length scales than temperature chaos. 
Right: Log-log plot of the same data, plotted as 1 — Q, to 
illustrate the quality of the scaling. 



have tested the temperature-dependence of the exponent 
and find that for T < 0.5 its value is practically inde- 
pendent of temperature, i.e., at T = 0.20 we are probing 
the low-temperature regime. 

Renormalization group arguments also suggest that 
the temperature and disorder chaos efi^ects are deeply 
related and characterized by the same universal scaling 
function 0, 0, [3] so that only nonuniversal pref- 

actors differ. In Fig. [21 we thus superimpose the data 
for both perturbations presented in Fig. [1] by rescaling 
ic- Using — 1.16 for both perturbations and multiply- 
ing ici^J) by a factor A « 8.7, we obtain a rather good 
superposition in the low-temperature region, and we con- 
clude that our data are thus compatible with the equality 



of the two scaling functions. From this -^-scale renormal- 
ization we also conclude that the length scale at which 
temperature chaos appears is approximatively ten times 
larger than the length scale needed for disorder chaos to 
appear, as has also been found in four space dimensions 
by Sasaki et al, see Ref. 0. Therefore temperature 
chaos is harder to probe than disorder chaos, as has al- 
ready discussed within a Migdal-Kadanoff approach on 
spin glasses (20i |. 

The superposition of the scaling functions shows devi- 
ations at larger temperatures. This is not surprising as 
assumptions made when deriving the scaling function are 
not valid for high T. For example, the Vt dependence 
of the entropy or the very existence of droplets is only 
valid for T < Tc [s^]. We thus believe that to perform 
a definitive test of universality of the scaling functions, 
very large system sizes at low temperatures with small 
temperature changes should be used. 

We also study the behavior of the scaling function. It 
decays as (Ic/L)'^/'^ for strong chaos (when £c/L < 1) [s^ . 
However, in the limit £c/L > 1, we obtain 1 — Q{x) oc 
x^^/^ which differs from the 1 — Q{x) cx x'^ behavior found 
at zero temperature in Ref. 32,]. This shows that the 
scaling function can have a more subtle behavior than 
what is naively expected from simple domain-wall argu- 
ments [H, [sl] . Note that the results remain unchanged 
if one takes the disorder average of the different overlaps 
in Eq. ^ independently, as done in Eq. (11) of Ref . [l3|. 

Finally, to better illustrate the mechanism of chaos, we 
study the distribution of the chaoticity parameter over 
the disorder for temperature chaos, i.e., we compute the 
chaoticity parameter Q as defined in Eq. ([4]) without the 
disorder average, and bin the data for different choices 
of the disorder to compute the distribution Pl{Qti,T2)- 
According to the droplet model, in the weak chaos regime 
(where ic > L), temperature chaos can manifest itself 
even on small length scales, but only for rare regions of 
space This means that even for small AT, when 

Q is very close to unity, the distribution is broad and 
rare samples with lower values are expected. Figure [3] 
shows the distribution Pt,{Qt,.t.,) for L = 10 for Ti = 
0.2 and different T2 = Ti + dJ: Even for modest 

AT, rare but large changes are clearly observed. This 
illustrates the weak chaos scenario presented in Ref. 17| : 
temperature chaos (at least in the weak regime) is not 
due to moderate changes in all samples, but rather due 
to larger changes in a few rare samples. 

Conclusions We have studied numerically disorder 
and temperature chaos in 3D Ising spin glasses and show 
that both disorder as well as temperature chaos are well 
described within a scaling/droplet description. In partic- 
ular, we find that the scaling variables have to be modi- 
fied as done in Eq. ([8|) when the difference in temperature 
is large. In addition, we show that the weak chaos regime 
is dominated by rare events where system-size droplets 
are fiipped. This has direct experimental implications be- 
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FIG. 3: (Color online) Distribution of the chaoticity param- 
eter over different realization of the disorder for L = 10 with 
Ti = 0.2 and r2 = Ti + AT 29]. In the case of small AT rare 
samples have very low values of Q while most of them remain 
unchanged ((3 = 1). The inset shows a linear- log plot of the 
data for the smallest values of AT for which Ti < T2 < 2Tc/3. 
For all T2 the distribution has a single peak, unlike for the 
random-energy models with entropic fluctuations [39l . |40| . 



cause the weak chaos regime has been argued to account 
in a quantitative way for the memory and rejuvenation 
effects Finally, we show that temperature and dis- 
order chaos might be described by similar scaling func- 
tions in the low-temperature regime, thus providing com- 
peUing evidence for the presence of a chaotic temperature 
dependence in spin glasses. This has also recently been 

431. This mechanism 



proven for mean-field systems [41 



is also responsible for step-size responses that could in 
principle be observed experimentally in mesoscopic sys- 
tems [4l|,|42|. Nevertheless, this behavior might change 



for larger system sizes and thus we propose to revisit 
the problem with better models [i^. Our findings will 
help interpret experiments on rejuvenation and memory 
effects in spin-glasses and other materials. 
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